Percolation transitions with nonlocal constraint 
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We investigate percolation transitions in a nonlocal network model numerically. In this model, 
each node has an exclusive partner and a link is forbidden between two nodes whose r-neighbors 
share any exclusive pair. The r-neighbor of a node x is defined as a set of at most N r neighbors of 
x, where N is the total number of nodes. The parameter r controls the strength of a nonlocal effect. 
The system is found to undergo a percolation transition belonging to the mean field universality 
class for r < 1/2. On the other hand, for r > 1/2, the system undergoes a peculiar phase transition 
from a non-percolating phase to a quasi-critical phase where the largest cluster size G scales as 
G ~ N a with a = 0.74 (1). In the marginal case with r = 1/2, the model displays a percolation 
transition that does not belong to the mean field universality class. 
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A networked system is in a percolating phase when a 
finite faction of nodes are interconnected via links to form 
a percolating giant cluster 0, Q ■ A percolating config- 
uration requires a macroscopic number of links, and the 
link density / given by the ratio of the number of links 
over the number of nodes is a control parameter for a 
percolation transition. It is well known that there is a 
threshold link density l c above which a percolating clus- 
ter appears. Recently, the percolation in complex net- 
works has been studied extensively for the understanding 
of emergent phenomena in complex systems [H, H| 

A prototypical model for percolation transitions is the 
random network model of Erdos-Rcnyi (ER) where each 
pair of nodes is connected independently and uniformly 
with the same probability Q. It displays a percolation 
transition at a critical link density l c = 1/2, and its criti- 
cal behavior exemplifies the mean-field (MF) universality 
class 0, Q . The universality class is extended by gener- 
alizing the model in various ways. For example, perco- 
lation transitions are studied in random networks with 
a power-law degree distribution 0, 0, @K in growing net- 
works [(I [l(J , in correlated networks [lj] - [l3| , and so on. 

In most network models, a linking probability between 
nodes depends only on a local information around in- 
volved nodes. Recently, a nonlocal model was proposed 
in Ref. (l4j . where a node is linked with a probability 
depending on the whole cluster size distribution. This 
model exhibits an intriguing explosive percolation tran- 
sition, and a lot of works have been performed to clarify 
its nature [l5l - [2l| . These studies suggest that nonlocality 
may be an important ingredient for the universality class 
of percolation critical phenomena. 

In this paper, we introduce a nonlocal network model 
under a so-called pair-exclusion (PE) constraint and in- 
vestigate its percolation transitions. Each node is as- 
signed to have its own exclusive counterpart. The PE 
constraint means that a link between two nodes is for- 
bidden if two clusters containing each node share any ex- 
clusive pair. Such a constraint was considered in evolving 
network models [22[ and mass aggregation models [23j |. 




FIG. 1: (Color online) Illustration of dynamics of the model 
with N = 16 and r = 1/2. Shaded regions represent the r- 
neighbors of nodes of species m, 713, and ns. One can add 
a link (represented with a dashed line) between nodes of m 
and 713, while a link (represented with a dotted line) between 
nodes of ni and 77,5 is forbidden because of an exclusive pair 
of species 712 represented with filled circles. 



In the context of evolving networks, the PE corresponds 
to a constraint against self- loops [22j . The PE type inter- 
action may also be relevant in a social group formation 
process among individuals with mutual conflicts [24| . We 
apply the PE constraint to a percolation problem. 

Consider a system consisting of ./V nodes. In order to 
impose the PE constraint, we assume that there are N/2 
particle species denoted by with i = 1, . . . ,N/2 with 
two members in each species (see Fig. [TJ. For each node 
x, a cluster C x is defined as a set of nodes that can be 
reached from x via links. We also define a r-neighbor 
V x as a subset of C x consisting of at most [N r ] nodes. 
Here < r < 1 and [N r ] denotes the integer part of N r . 
It includes nodes in the ascending order of the distance 
from x. If the size of C x is less than [N r ], then V x is 



2 




simply given by C x . So, the sizes of V x and C x denoted 
by \V X \ and \C X \, respectively, satisfy 

\V x \=xmn{\C x \,[N r ]} . (1) 

Initially, one starts with N isolated nodes with no link. 
Then, links are added sequentially in the following way: 
One selects two nodes x and y at random. Two nodes 
are connected with a link if their r-neighbors V x and 
V y do not share any exclusive pair. Otherwise, linking 
is rejected (see Fig. [1]). The parameter r represents a 
strength of the nonlocal constraint. 

We are interested in how the onset of a percolation 
and the scaling of a percolating cluster are affected by 
the nonlocal constraint parameterized by r. The relevant 
measure is the sample-averaged largest-cluster size G as 
a function of the number of added links L or the link 
density I = L/N. The percolation order parameter is 
given by o = G/N. The mean cluster size S is also 
measured Q- 

Figure [2] shows the overall behavior of the order pa- 
rameter measured in systems with N = 64000 nodes and 
averaged over 10000 samples. The order parameter dis- 
plays a sharp transition at small values of r, while it 
displays a weak transition behavior at large values of r. 
Apparently the nonlocal constraint suppresses the emer- 
gence of a percolating cluster. Percolation properties in 
both cases will be studied in detail. 

When r = 0, the PE constraint becomes purely lo- 
cal in that a linking trial is rejected only when selected 
nodes are an exclusive pair. There are N/2 exclusive 
pairs among N(N — l)/2 possible selections. So a linking 
is rejected with the probability 

p = l/(N-l), (2) 

which can be ignored in the N — > oo limit. Hence, the 
system becomes equivalent to the ER random network 
model in the N — > oo limit, and undergoes a percolation 
transition in the MF universality class at I = l c = 1/2. 



With nonzero r, a linking probability P xy for ran- 
domly selected nodes x and y depends on the size of 
their r-neighbors V x and V y . There are \V X \ ■ \V y \ pos- 
sible combinations between nodes and a linking trial is 
accepted only when none of them are an exclusive pair. 
Hence, the linking probability is given by 

P X ,y = (l-P0) lVMl (3) 

with po in Eq. It is approximated as 

P Xty ^ e -WMV y \/N (4) 

in the large N limit. 

Equation j4| gives a hint on the role of the nonlo- 
cal constraint. We note that jV^I < \C X \ from Eq. |T]). 
Hence, nodes belonging to a larger cluster have a higher 
rejection probability. Consequently, the nonlocal con- 
straint suppresses the growth of large clusters. We also 
note that \V X \ < [N r ] from Eq. (TT|). So, the quantity in 
the exponent of Eq. ^ is bounded as 

\V x \-\V y \/N<N^ (5) 

for any pair of x and y. It indicates that the nonlocal con- 
straint leads to a different effect depending on whether 
r < 1/2 or not. 

Firstly, we consider the case with r < 1/2. It will be 
referred to as a weak constraint region. In this region, 
the linking probability in Eq. ^ converges to 1 with a 
small correction of the order of 0(7V~( 1 ~ 2r )) at most. 
Therefore, the nonlocal constraint is negligible and the 
model is expected to display the percolation transition 
belonging to the MF universality class. Furthermore, the 
percolation threshold is equal to that of the ER random 
networks, l c (r) = 1/2. 

In order to characterize the percolation transition, we 
have measured numerically the order parameter g, the 
density of the largest cluster, and the mean cluster size 
S. They satisfy the finite-size-scaling (FSS) forms 

g(l,N) - N-^F g {{l-l c )N^) , (6) 
S(l,N) = N^-TsHl-QN 1 /") (7) 
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with the order parameter exponent /3, the mean clus- 
ter size exponent 7, and the FSS exponent v 0. The 
MF universality class is characterized with the exponents 
Pmf = 1, Imf = 1, and vmf = 3 0. Figures E£a) and 
(b) present the scaling plots for g and S at r = 0.3. Each 
data set collapses well onto a single curve with l c = 1/2 
and the MF critical exponents. We obtained the same 
scaling behavior for other values of r < 1/2. This con- 
firms our expectation based on the scaling of the linking 
probability. 

Secondly, we consider the case with r > 1/2, which will 
be referred to as a strong constraint region. Here, the 
scaling behavior of P x _ y in Eq. ^ depends on the cluster 
size distribution. Since the PE constraint suppresses the 
growth of the largest cluster, its size G is bounded above 
by the largest cluster size of the ER random network. To 
a given value oil < 1/2, the ER random network displays 
a subcritical scaling G ~ In TV Q. Consequently, the 
linking probability P xy for any x and y converges to 1 
with a correction of the order of 0((ln TV) 2 /TV) at most. 
Therefore, for I < 1/2, the system becomes equivalent 
to the ER random network model and lies in the non- 
percolating phase with g ~ In N/N. 

In the strong constraint case, the constraint is irrel- 
evant for I < 1/2 since \V X \ < \C X \ < G = In TV for 
all x. It becomes relevant when the largest cluster size 
reaches a scaling G ~ TV 1 / 2 . One can locate the thresh- 
old value l c {r) at which G follows a critical power-law 
scaling from the following scaling argument: The sub- 
critical scaling G ~ In TV for I < 1/2 in the previous 
paragraph implies that l c {r) > 1/2. In the ER random 
network model (r = 0), the largest cluster follows the 
critical scaling G ~ n^-Pmf/^mf _ /y2/3 at itg perco i a _ 

tion threshold I = 1/2 Q. The exponent value is larger 
than 1/2. This implies that the nonlocal constraint is 
already relevant at I = 1/2 suggesting that l c (r) < 1/2. 
Therefore, the threshold value is given by l c (r) = 1/2 at 
all values of r. Furthermore, we expect that the critical 
largest cluster size follows a power-law scaling 



G C (TV) - TV Q <= 



(8) 



with l/2<a c <2/3in the strong constraint case. 

We have measured numerically the largest cluster size 
G at I — l c to estimate the exponent a c . Figure SJa) 
shows the numerical data for the effective exponent 
a c (TV) ee ln(G c (2TV)/G c (TV))/ln2 at several values of r. 
We find that 



0.56(5) 



(9) 



for r > 1/2. The exponent value is different from the MF 
value (1 — Pmf/vmf) = 2/3. Note that a c ~ 2/3 in the 
weak constraint region. The marginal case with r =1/2 
will be considered later. 

Can the system have a macroscopic percolating cluster 
of size G = 0(N) in the strong constraint case (r > 1/2)? 
Let us recall the mechanism leading to a percolating clus- 
ter in the random network model (r = 0). As one adds 
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FIG. 4: (Color online) (a) Effective exponent a c (N) for the 
critical power-law scaling of G c ~ TV Q ° at the percolation 
threshold I = 1/2 at several values of r as indicated by leg- 
ends, (b) Effective exponent a (TV) for the power-law scaling 
of G ~ TV" at r = 0.3 (filled symbols) and r = 0.5 (open sym- 
bols), (c) Effective exponent a (TV) for the power-law scaling 
of G ~ TV" at r = 0.7 (filled symbols) and r = 1.0 (open 
symbols). In (b) and (c), link densities are I — 0.75 (o), 
I = 1.0 (□), I = 1.25 (o), and I = 2.0 (A). 



links between pairs of nodes, clusters merge with each 
other and grow. Nodes are selected for linking at ran- 
dom with the same probability. However, growth rates 
of clusters are not uniform. Larger clusters grow faster 
than smaller ones because a cluster is chosen for linking 
with the probability proportional to its size. Hence, once 
there emerges a dominant cluster whose size is consid- 
erably larger than the others, it grows even faster and 
eventually forms a giant cluster. 

The strong constraint suppresses the growth of large 
clusters severely. Especially, the linking probability be- 
tween clusters of size 0(N q ) with q > 1/2 is vanishingly 
small. It suggests that the strong constraint does not 
allow for a single dominant cluster. Then, a possible sce- 
nario is that there appears a set of mesoscopic clusters 
of a same FSS behavior. It is reminiscent of a so-called 
powder keg scenario for the explosive percolation [2o| . 
Explosive percolation models also have a nonlocal con- 
straint suppressing the growth of large clusters. It results 
in a powder keg of large clusters, which suddenly merge 
into a giant cluster at a percolation threshold. The con- 
straint in our model is much stronger. For example, once 
a linking between two nodes is rejected, they are never 
connected afterwards. It suggests that a giant cluster of 
size G = 0(N) may be improbable. Detailed properties 
are investigated numerically. 

The size of the largest cluster G at I > l c is measured 
numerically. It turns out to follow a power-law scaling 



G - TV" 



(10) 



Numerical values of the exponent are estimated by us- 
ing an effective exponent a(N) = ln[G(2TV)/G(TV)]/ln2. 
Figures 0|b) shows the plot of a(TV) at r = 0.3. It con- 
verges to 1 for all values of I > l c indicative of a macro- 
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FIG. 5: Phase diagram 



scopic percolating cluster. On the other hand, Fig. @|c) 
shows that the exponent converges to 



a = 0.75(1) 



(11) 



at all values of r > 1/2. Therefore we conclude that the 
strong constraint leads to a quasi-critical phase charac- 
terized with the exponent a ~ 0.75. 

In the marginal case with r = 1/2, the system displays 
a distinct percolation transition. The largest cluster size 
scales as G ~ N a <= with a c = 0.61(2) (see Fig.gja)) at 
the percolation threshold and there exists a macroscopic 
giant cluster with G ~ N 1 in the supercritical phase (see 
Fig.SJb)). Using the FSS analysis based on Eqs. (j6|) and 
0, we find that the transition is characterized with the 



MF critical exponents (see Figs.[3jc) and (d)) 



0/D = 0.38(2), -y/D = 0.28(3), 1/V = 0.31(3) . (12) 

They satisfy the scaling relation 2f3/D+j/D = 1 Q within 
error bars. These non-MF exponents indicates that the 
marginal case constitutes a distinct universality class for 
percolation transitions. 

In summary, we have introduced a percolation model 
with a nonlocal constraint and investigated the nature 
of percolation transitions numerically. Our results are 
summarized in the phase diagram in Fig. [5l When 
I < l c = 1/2, the system is always in the non-percolating 
phase where G ~ In N. When I > l c , the system is in the 
percolating phase with G ~ TV 1 for r < 1/2 while it is in 
the quasi-critical phase with G ~ N a with a ~ 0.75 for 
r > 1/2. The critical line I = l c in the weak constraint re- 
gion (solid line) belongs to the MF universality class. The 
critical line I — l c in the strong constraint region (dotted 
line) separates the non-percolating phase and the quasi- 
critical phase. The marginal case at r = 1/2 belongs to a 
distinct universality class characterized with the non-MF 
exponents in Eq. (fT2j) . Our study shows that the nonlocal 
constraint leads to rich percolation critical phenomena. 
It calls for further study of percolation transitions in sys- 
tems with nonlocal constraints. 
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